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We present measurements of the angular cross-correlation between luminous red galaxies from 
the Sloan Digital Sky Survey and the cosmic microwave background temperature maps from the 
Wilkinson Microwave Anisotropy Probe. We find a statistically significant achromatic positive 
correlation between these two data sets, which is consistent with the expected signal from the late 
Integrated Sachs- Wolfe (ISW) effect. We do not detect any anti-correlation on small angular scales 
as would be produced from a large Sunyaev-Zel'dovich (SZ) effect, although we do see evidence for 
some SZ effect for our highest redshift samples. Assuming a flat universe, our preliminary detection 
of the ISW effect provides independent physical evidence for the existence of dark energy. 

PACS numbers: 98.65.Dx,98.62.Py,98.70.Vc,98.80.Es 



I. INTRODUCTION 

As photons from the Cosmic Microwave Background 
(CMB) travel to us from the surface of last scatter- 
ing, they can experience a number of physical processes. 
These include the Sunyaev-Zel'dovich (SZjS^I) effect, 
which is the inverse Compton scattering of CMB photons 
by a hot (T > 10^ K) ionized_gas0, and the late Inte- 
grated Sachs-Wolfe (ISW^lllll) effect, which is the in- 
tegrated differential gravitational redshift caused by the 
evolution of gravitational potentials along the path trav- 
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eled by the photons. At frequencies less than 217 GHz, 
the SZ effect produces an anti-correlation between the 
temperature of the CMB and galaxies on small angular 
scales. The ISW effect generates a weak positive cor- 
relation at large angular scales for all freauenciesp^. 
For a universe with a flat geometry [23, the detection 
of any ISW effect provides strong physical evidence for 
the existence of dark energy, and can be used to mea- 
sure its equation of state [a, IS 0|- In this paper, 
we attempt the detection of the ISW and SZ effects 
through the cross-correlation of new, large area, clean 
galaxy and CMB maps of the sky 9, 22], which reduce 
the errors due to sample variance on large scales and the 
Poisson error on small scales relative to that of previous 
attempts H in USES- 
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The latter is a new, hierarchical equal-area pixelization 
scheme developed specifically for the SDSS geometry 
and fast correlation statistics. For our primary calcula- 
tions we use the higher resolution SDSSPix for accurate 
weighting, while for the extensive Monte Carlo simula- 
tions discussed below, we use HEALPix to take advan- 
tage of the SpICE'2^ for fast cross-correlation estima- 
tions. Both methods required re-pixelization of the orig- 
inal maps. We have tested our re-pixelization methods 
by comparing the auto-correlation functions of the orig- 
inal and re-sampled maps and we find identical results 
and errors. 

For the cross-correlation measurement, we use a simple 
pixel-based estimator. For each pixel i in our map, we 
determine the galaxy and temperature over-densities (A^ 
& A"^, respectively). 



FIG. 1: The photometric redshift distributions for the four 
LRG subsamples, taking into account the covariance between 
photometric redshift and galaxy type. To select these galax- 
ies, we imposed a color cut of 0.7{g — r) + 1.2((r — i) — 0.18) > 
1.6 and (g — r) > 1. We then define dx = (r— «) — ((/— r)/8, and 
step through d± in increments of 0.2 starting at 0.2 < d± < 
0.4 and ending at dx > 0.8 to determine our four redshift 
samples, from lowest to highest redshift respectively. There 
are 0.4 million LRGs in the z ~ 0.35 subsample, 0.8 million in 
the z ~ 0.43 subsample, 1 million in the z ~ 0.49 subsample, 
and 0.7 million in the z ~ 0.55 subsample. 



II. DATA & ESTIMATORS 

The parent SDSS galaxy data used herein is nearly 
25 million galaxies taken from ~ 3400 square degrees 
of the sky i mag ed by the Sloan Digital Sky Survey 
(SDSS 0, 111 111 111 111 111 111 El)- We apply a flux 
limit of i < 21 and regions with poor image quality have 
been excluded from our analysis {i.e., data with a point- 
spread function of worse than 1.5 arcseconds). From this 
basic data set, we have constructed four different subsets 
of galaxies using both the Luminous Red Galaxy (LRG) 
selection methodjlO| and a photometric redshift estimate 
for every galaxyj^- As shown in Figure^ these subsam- 
ples span a range of redshift from z ~ 0.3 to z ~ 0.8, 
providing an efficient method for tracing the large-scale 
distribution of galaxies in the Universe as a function of 
redshift. 

For the CMB data, we use the 3 primary combined- 
channel maps (Q, V, W) from the Wilkinson Microwave 
Anisotropy Probe (WMAP'27'l), and the "clean" map of 
Tegmark et al. ^ . We also apply a mask to exclude fore- 
ground contamination from the Galaxy (kpl2 for the V 
& W bands and kp2 for the Q band) . To isolate the IS W 
effect, we also generate a smoothed version of the "clean" 
map, by applying a Gaussian kernel with a FWHM of one 
degree to remove CMB primary anisotropics on small an- 
gles. 

The natural pixelization scheme for the WMAP maps 
is HEALPix while for the galaxies it is SDSSPix|; 



Af = 



T,-T, 



(1) 



where Ti is the temperature in pixel i, T is the mean 
temperature, is the total number of galaxies in pixel 
I, and nSj is expected number of galaxies in pixel i. In 
addition, we can optimally weight each pair of pixels to 
account for the variation of signal-to- noise (S/N) in the 
pixels. Therefore, the two-point cross-correlation func- 
tion at angular separation 9^ is given by 



(2) 



where 9a is the central value of an angular correlation bin, 
fi is the fraction of pixel i unmasked, is the inverse of 
the noise in the WMAP pixel, and Qfj is unity if the 
separation between the pixels is within this angular bin 
and zero otherwise. On large angular scales the optimal 
CMB weighting is uniform (i.e., st = 1, except where 
masked) as the signal is dominated by sample variance 
rather than Poisson noise. For completeness, we calculate 
the cross-correlation functions for both the uniform and 
noise-weighted schemes. 

To estimate the covariance C{9a , 9f}), we use two meth- 
ods: jack-knife errors and random CMB maps. For the 
former, we divide our total area into Nj equal-area sub- 
samples and use the jack-knife estimator ^25.J. 
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where Wgf{0) is the mean of WgT:{9) for the Nj mea- 
surements and WgT.i(0) is the measurement of the cross- 
correlation excluding the zth subsample. In order to pro- 
duce a stable estimation of the covariance matrix (non- 
singular and positive definite), we set Nj ~ 2Ne, where 
Ne is the number of angular bins. Tests with different 
values of Nj near to 2Nq generally produced values of 
(see Section FlV B|l within 15% of the values obtained 
using Nj = 2Ng. 



While the jack-knife method gives a good approxima- 
tion of the sample variance, it is less sensitive to variance 
on angular scales larger than the survey size. For a more 
conservative estimate of the covariance, we generated 200 
random CMB maps based on the angular power spec- 
trum of the "clean" map and 200 random CMB maps 
based on the angular power spectrum of our smoothed 
WMAP map. These maps were cross-correlated against 
the four galaxy maps (using SpICE0|) and the results 
were re-sampled into the same angular bins as used in 
the SDSSPix analysis. The covariance matrix is then 



7V=200 



i=l 



x{w^Ti9l3) - WgT,i(6'/3)). 



(4) 



III. RESULTS 



We present the cross-correlation functions for the four 
redshift subsamples and the four WMAP maps in Fig- 
ure |2| For the three highest redshift galaxy subsamples, 
we find a consistent non-zero signal. We also find that 
the cross-correlation signal is achromatic in the Q, V, 
W frequency maps. We see little evidence for an anti- 
correlation between the galaxy density and CMB tem- 
perature on the smallest angular scales, as expected from 
the SZ effect. For the lowest redshift subsample, we see 
a different shape to the cross-correlation function and 
all the angular bins are consistent with zero within their 
errors. This lowest redshift subsample also shows a lack 
of agreement between the Q, V, and W frequency maps. 



IV. DISCUSSION 
A. Systematic Error Tests 

Our primary concern is the possible contamination 
of the CMB signal by residual Galactic synchrotron or 
dust emission and the galaxy catalogs by Galactic stars. 
This could result in a non-zero correlation between the 
maps. To check against this, we calculate the stellar over- 
density in each pixel, using the same color cuts as for 
the four galaxy subsamples discussed above, and substi- 
tute these over-densities into Equation |21 For the three 
highest redshift subsamples, we find flat, featureless star- 
WMAP cross-correlation functions consistent with zero 
at the la level. For the lowest redshift sample (z ^ 0.35), 
the star- WMAP cross-correlation function is nearly iden- 
tical to the galaxy- WMAP measurement, suggesting sig- 
nificant contamination. 

The cross-correlation functions for the three high- 
est redshift galaxy subsamples are achromatic in the 
WMAP Q, V and W frequency maps. In addition to 
the visual agreement in Figure [3 we cross-correlated 
the three highest redshift galaxy subsamples against four 



2 
1.5 
1 

0.5 


.-0.5 

J) 

1.5 
1 

0.5 


-0.5 



T — I I I I I l| 1 1 — I I I I ll| 

z ~ 0.57 



T — I I I I I l| 1 1 — I I I I 1 1| 

z ~ 0.49 




I Lj ■ 



9 (degrees) 



FIG. 2: The uniform- weighted galaxy- WMAP cross- 
correlation functions for the four galaxy subsamples. The 
shaded boxes show the measurement and the la jack-knife 
errors for the "clean" galaxy-WMAP cross-correlation func- 
tion. The errors on the other functions are nearly identical. 
The Q map is given by the red line, V by the green line, W 
by the blue line and smoothed map by the cyan line. 



linear combinations of the Q, V, and W maps. These 
combination were designed to detect achromatic effects 
(Q — V and V — W) and possible contamination from 
synchrotron ((9Q — 12V -|- 3W)/16) and dust emission 
((3Q - 12V + 9W)/16). In all 12 cross-correlation func- 
tions, the signal was suppressed relative to our original 
galaxy-WMAP signal by a factor of ~ 10 on all angular 
scales with errors consistent with Poisson noise from the 
galaxies. Likewise, the shapes of the cross-correlation 
functions were dissimilar to those seen for the galaxy- 
WMAP cross-correlation functions. These observations 
confirm that our galaxy-WMAP measurements are not 
due to contamination by synchrotron or dust emission 
from the Galaxy, or synchrotron emission the host galax- 
ies themselves |2(](. i.e., over the frequency range probed 
by WMAP, the intensity of any synchrotron emission 
should decrease by an order of magnitude from Q to W. 



B. Significance Tests 

As with any angular correlation measurement, the in- 
dividual angular bins are highly correlated. Therefore, 
we need to consider the the full covariance matrix (Cn^j) 
when checking the statistical significance of our mea- 
surements against the null hypothesis of a zero cross- 
correlation function. Our values are given by 



(5) 
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TABLE I: The values derived from both covariance ma- 
trices (jack-knife errors given in parentheses) for A'^e — 10 
angular bins 



Subsample 


Q 


V 


W 


z ~ 0.55 
z ~ 0.49 
z ~ 0.43 
z ~ 0.35 


10.4 (16.3) 
19.9 (28.3) 
8.4 (16.3) 
33.9 (14.6) 


18.6 (20.4) 
21.6 (32.0) 
5.9 (19.5) 
54.2 (24.6) 


35.7 (25.3) 
19.9 (17.0) 
13.0 (37.2) 
27.3 (33.3) 


Subsample 


"Clean" 


Smoothed 




z ~ 0.55 
z ~ 0.49 
2 ~ 0.43 
2 ~ 0.35 


15.8 (27.1) 
17.0 (30.2) 
6.8 (34.1) 
32.5 (27.3) 


7.0 (20.8) 
11.0 (24.8) 
5.6 (30.1) 
12.8 (11.2) 





where Cr^j is the covariance matrix derived from either 
Eqn^(random maps) or Eqn|21 (jack-knife errors) respec- 
tively. The values for each of our uniform weighted 
cross-correlation functions, including the smoothed map, 
are given in Table ^ using both covariance matrices. 
These values are in reasonable agreement with other mea- 
surements of similar cross-correlation functions at larger 
angles 0,1^. The choice of weighting for the CMB pixels 
does not change our results significantly. 

For the three highest redshift samples, we exclude the 
null hypothesis at > 90% confidence for every WMAP 
band and the "clean" and smoothed maps using the jack- 
knife errors. For the random map errors, we can exclude 
the null hypothesis at > 90% for only the two highest 
redshift samples (excluding the Q band and smoothed 
map). Excluding angular bins on scales < 2 degrees in- 
creases the significance of the smoothed detection with 
the random map errors but not to the 90% confidence 
threshold. Checking the number of random maps with 
greater than the values in Table ^ verifies that our 
confidence estimates are correct. 

We complement the statistical a naly sis by apply- 
ing the False Discovery Rate (FDR [la IS) technique to 
the combination of all of our cross-correlation functions 
into a single significance test of our whole signal. The 
utility of the FDR statistic over combining the individ- 
ual x^'s for different redshift slices comes from the fact 
that it is conservative and designed to work with highly 
correlated datasets. FDR works by controlling a, the 
fraction of false detections compared to the total num- 
ber of detections^. For a given a, the FDR technique 
provides a threshold above which bins are rejected from 
a null hypothesis, in this case Wgx^i = for all i angular 
bins, with the guarantee that the true number of errors is 
always < a. Using the three highest redshift subsamples, 
cross-correlated with the five CMB maps (5x10x3= 150 
bins in total), we reject, with an a = 0.25, 144 bins with 
the jack-knife errors and 113 bins for the random CMB 
map errors. Therefore, at least 108 (84) of our 150 angu- 
lar bins are rejected from the null hypothesis using the 
jack-knife (random maps) errors. 



V. THEORETICAL MODELS 

We now compare our results to simple physical models 
of the ISW and SZ effects. A more detailed theoretical 
analysis of our results will appear in a future paper. We 
assume here a ACDM cosmology given by the fits to the 
WMAPext data[23: ^ot = 1, /i = 0.72, n^, = 0.29, 
fib = 0.047, 0-8 = 0.9, and = 0.99. 

We expect our observed cross-correlation functions to 
be the sum of the induced correlations from the ISW and 
SZ effectsii, i.e., 



WgriO) = wisw(^) + wsziO). 



(6) 



The ISW effect is generated by the decay of gravitational 
potentials in an open or dark energy dominated universe. 
This decay couples to the dark matter momentum, mean- 
ing that our galaxy-temperature cross-correlations will be 
a function of the galaxy-momentum power spectrum: 



wiswiO) = b 



dh 

y / dxK.isw{k,e,x)D{x)P{k), (7) 



where x is the co-moving distance, D{x) is the linear 
growth factor, P(k) is the linear dark matter power spec- 
tra and h is the linear galaxy bias at the peak of the 
galaxy redshift distribution (z). The kernel (/Cisw) is 
given by 



D{x) 



a(x) 



Mkex), 



(8) 

where Hq is the Hubble constant, Tq is the CMB tem- 
perature today in /iK, Wg is the normalized galaxy dis- 
tribution and Jo is the zero order Bessel function. 

The SZ effect is the cross-correlation of the average 
free-electron pressure (oc nJTc) along the line of sight, 
with the projected galaxy density. The pressure-galaxy 
power spectrum is therefore. 



wsz{e) ^ n,T,hpb I dkk I dxlCsz{k,6,x)D\x)P{k), 

(9) 

where n^, To, and 6p arc the electron density, electron 
temperature, and pressure bias of the gas respectively, 
while represents the path length weighted average. In 
the Rayleigh- Jeans portion of the thermal CMB spec- 
trum probed by WMAP, the kernel (A^sz) is given by 



/Csz(fc,0,x) = -2 



o-TfcBToWg(x) 



MkOx), (10) 



where ctt is the Thompson scattering cross-section and 
fee is Boltzmann's constant. If we set rio to the mean 
electron density at z, we have two free parameters to fit 
to the data, i.e., b and Tcbp. 

Figure |3| presents an example of our best fit mod- 
els for the z ^ 0.49 galaxy subsample cross-correlated 
with the W-channel WMAP map. For this cross- 
correlation function, we find that the best fit model is 
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FIG. 3: The comparison of our theoretical predictions to our 
measurement of the W-channel WMAP map cross-correlated 
with the z ~ 0.49 galaxy subsample. The fitted ISW cross- 
correlation function is given by the red line, the SZ by the 
green line and the sum of the two is given by the blue line. 
We find b = 5.47 ± 1.82 and TJ^ = 0.15 ± 0.61, where the 
errors are unmarginalized. 



VI. CONCLUSIONS 

We present here measurements of the cross-correlation 
function between the currently available SDSS galaxy 
data and the WMAP CMB maps of the sky. The in- 
dividual measurements in Figure El are significant at the 
> 90% confidence, while FDR strongly rejects the null 
hypothesis. The achromatic nature of the signal (in the 
Q, V, and W frequency channels) confirms the signal is 
of cosmological origin. Our results are consistent with 
the NVSS-CMB correlations |1 l2l|. but may be harder 
to reconcile with other optical-CMB measurements of 
the ISW effect [nl 113. These results are consistent with 
the predicted ISW effect (which is largely independent 
of dark energy properties) and a minor SZ contribution. 
Assuming a flat universe, our detection of the ISW effect 
provides independent physical evidence for the existence 
of dark energy. In future papers, analyzing all available 
SDSS data, we will present a more detailed compari- 
son between the theoretical models and the data, as well 
consistency checks of the observed correlation amplitude 
with galaxy biases and dark energy models. 



preferred to the null hypothesis (Section IIVB|I at the 
99% confidence level, based on the observed difference 
between the values of the model fit (Xmodei) ^^'^ 
null hypothesis ixluu): *-e-> = xlnodei - xluU ^ 
9.1 for 2 degrees of freedom, using the jack-knife er- 
rors. If we consider the three highest redshift galaxy 
subsamples, cross-correlated with the 5 WMAP maps 
(W,Q,V," clean" , smoothed), we find that, using the jack- 
knife covariance matrix (Cj), 14 of these 15 cross- 
correlation functions are better fit by our physical model 
than by the null hypothesis at the > 95% confidence level 
(based on the Sy^ values). If we assume the more conser- 
vative random map covariance matrix (C^), then four of 
the fifteen functions are again rejected at the same con- 
fidence (or eight at > 80% confidence). In summary, a 
majority of our cross-correlation functions are much bet- 
ter fit by our physical model than by the null hypothesis. 

We find a somewhat stronger SZ signal in the higher 
redshift bins, although the SZ effect is always sub- 
dominant to the ISW signal. We also see higher values 
for h for the z ^ 0.55 galaxy subsample than for z 0.43 
galaxy subsample, but the fit values of b are sensitive to 
the details of our photometric redshift distributions. We 
will present a more detailed analysis and interpretation 
in a future paper. 
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